Support for the thermal origin of the Pioneer anomaly 



(N 

O 

(N 

< 



Slava G. Turyshev\ Viktor T. Toth^, Gary Kinsella\ Siu-Chun Lee^ Shing M. Lok^ and Jordan Ellis^ 

^ Jet Propulsion Laboratory, California Institute of Technology, 
4800 Oak Grove Drive, Pasadena, CA 91109-0899, USA 

'^Ottawa, ON KIN 9H5, Canada and 

^Applied Sciences Laboratory, 13111 Brooks Drive, Suite A, Baldwin Park, CA 91706-7902, USA 

We investigate the possibility that the anomalous acceleration of the Pioneer 10 and 11 space- 
craft is due to the recoil force associated with an anisotropic emission of thermal radiation off the 
vehicles. To this end, relying on the project and spacecraft design documentation, we constructed 
a comprehensive finite-element thermal model of the two spacecraft. Then, we numerically solve 
thermal conduction and radiation equations using the actual flight telemetry as boundary condi- 
tions. We use the results of this model to evaluate the effect of the thermal recoil force on the 
Pioneer 10 spacecraft at various heliocentric distances. We found that the magnitude, temporal 
behavior, and direction of the resulting thermal acceleration are all similar to the properties of the 
observed anomaly. As a novel element of our investigation, we develop a parameterized model for 
the thermal recoil force and estimate the coefficients of this model independently from navigational 
Doppler data. We find no statistically significant difference between the two estimates and conclude 
that once the thermal recoil force is properly accounted for, no anomalous acceleration remains. 

PACS numbers: 04.80.-y, 95.10.Eg, 95.55.Pe 
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INTRODUCTION 

The anomalous acceleration of the Pioneer 10 and 11 
spacecraft [l[ is a discrepancy between modeled and ob- 
served radio-metric Doppler data received from the two 
vehicles. The discrepancy can be eliminated by incorpo- 
rating a constant sunward acceleration of unknown origin 
with magnitude of ap = (8.74 ± 1.33) x 10"^° m/s^ into 
the orbital model 0] . The presence of this acceleration, 
which became known as the Pioneer anomaly, was seen 
by many as a breakdown in the gravitational inverse- 
square law of general relativity that reveals itself in the 
dynamics of the outer solar system. There were also pro- 
posals that this anomaly was due, at least in part, to the 
waste heat emitted by the spacecrafts' radioisotope ther- 
moelectric generators (RTGs) , their electrical subsys- 
tems 3 , or both , motivating a thorough investigation 
of the spacecraft systematics (for review, see 0). 

Our current investigation began with the recovery of 
the entire telemetry record of both spacecraft, substan- 
tial project documentation, and additional Doppler data 
The analysis of the extended Doppler data set was 
completed recently Q, confirming the existence of the 
anomaly. This analysis also showed that although the 
direction of the anomalous acceleration cannot be deter- 
mined unambiguously, the Doppler residuals improve if 
one considers a temporally varying acceleration in the 
direction of the Earth, which is consistent with earlier, 
similar results 0,[l3|- These results are suggestive: If the 
acceleration were due to thermal recoil force, it would be 
along the spacecraft spin axis, which generally points in 
the direction of the Earth and would have temporally 
decreasing magnitude consistent with the decay of the 
on-board radioisotope fuel (Pu^'^®). 



In this Letter we report on the completion of a finite- 
element thermal model of the Pioneer 10 spacecraft. This 
model relies extensively on the project and spacecraft 
design documentation and was validated using redun- 
dant flight telemetry data. A parameterized form of the 



thermal force model |11| was also incorporated into a 
Doppler analysis, which yielded an independent estimate 
of coefficients characterizing the thermal recoil force. We 
compare the outputs of the two independent analyses - 
Doppler and thermal - to show that no statistically sig- 
nificant anomalous acceleration remains in the data. 



FINITE-ELEMENT THERMAL MODEL 

A comprehensive finite-element (FE) thermal model 
(Fig. [ij of the Pioneer 10 and 11 spacecraft was con- 
structed at the Jet Propulsion Laboratory (JPL) in col- 
laboration with the Applied Sciences Laboratory (ASL). 
The geometric and thermal models of the spacecraft were 
constructed using the SINDA/3D thermal modeling soft- 
ware [HI. While the software provides the capability 
to build a numerical model directly from CAD drawing 
files, no such files exist for a spacecraft designed 40 years 
ago. For this reason, the model was built in a more te- 
dious manner by specifying the coordinates of the ver- 
tices of each modeled spacecraft surface, using available 
blueprints and recovered project documentation. The 
spacecraft geometric model was built with a Monte Carlo 
based radiation analyzer (TSS) to calculate the radia- 
tive exchange factors using infrared emittance values for 
modeled surfaces specified within it. The model incorpo- 
rated approximately 3,300 surface elements, 3,700 nodes, 
and 8,700 linear conductors. The spacecraft thermal- 
mechanical configuration is simulated by a network of 
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FIG. 1: Illustrative representation of the thermal model of the 
Pioneer 10 spacecraft evaluated at 40 AU. Top left: spacecraft 
body interior (temperature range: blue —16° C, red +10° C); 
Bottom left: spacecraft exterior (blue —155° C, red —108° C); 
Right: entire spacecraft (blue -213° C, red +136° C). Un- 
modeled struts that connect the RTGs to the spacecraft body 
are indicated with yellow-black dashed lines. 

TABLE I: Pioneer 10 telemetered power at select heliocentric 
distances. Externally vs. internally located components are 
indicated where applicable (only 5 out of 11 distances shown). 



Power (W) 



Description 


3 AU 


10 AU 


25 AU 


40 AU 


70 AU 


Science, internal 


12.6 


12.6 


11.9 


8.8 


0.8 


Science, external 


6.4 


8.4 


6.4 


6.4 


0.0 


Subsystems 


20.2 


20.2 


20.4 


20.2 


19.5 


Electrical, internal 


63.2 


46.3 


35.5 


28.4 


17.5 


Electrical, external 


8.1 


4.7 


2.7 


2.3 


0.1 


TWT" thermal 


18.6 


18.6 


19.8 


19.5 


21.2 


Transmitter 


9.2 


9.2 


8.0 


8.3 


6.6 


Total consumed 


138.3 


120.0 


104.7 


93.9 


65.8 


RTG generated 


148.5 


127.1 


107.1 


94.0 


67.2 


Cable loss 


6.9 


5.3 


4.0 


3.2 


1.7 


Total available 


141.6 


121.8 


103.1 


90.8 


65.5 


Difference 


+3.4 


+1.8 


-1.6 


-3.1 


-0.4 



"Traveling wave tube amplifier. 



thermal capacitance, conductive couplings, and radiative 
exchange factors between all surfaces and to deep space. 

The software numerically solves the energy equation 
using equipment power dissipation from the spacecraft 
fight telemetry records (see TablelJ). RTG power was esti- 
mated using the well-known half-life, r = 87.74 yr, of the 
238pu radioisotope fueh Q^tgW = [2-^'-'"^/^]Q,tg{h), 
where io = July 1, 1972 and Qrtg(io) = 2578.179 W. The 
objective was to calculate the temperature distribution of 
all spacecraft surfaces. To accommodate the limitations 
imposed on us by thermal modeling software, the angular 
distribution of the radiative emission from the spacecraft 
to space was calculated by solving the thermal radiation 
equations with the spacecraft positioned at the center of 
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FIG. 2: Pioneer 10 thermal power contributing to sunward 
acceleration, including solar heating and reflected solar ra- 
diation (hollow circles), solar heating only (hollow squares) 
with reflected solar radiation removed, and all solar effects 
removed (filled circles). Positive values indicate radiation di- 
rected away from the Sun, resulting in a sunward acceleration 
of the spacecraft. 

a large (i.e., with radius of 40 high-gain antenna, or HGA, 
diameters), black spherical control surface. The amount 
of spacecraft radiative emission absorbed by each element 
of this control surface corresponds to the amount of mo- 
mentum carried in this direction. 

The model incorporated some parameters the values of 
which were less well known: e.g., the effective emissivities 
of multilayer insulation blankets or conductive couplings 
of certain structural elements. Redundancies in the flight 
telemetry were used to refine the estimates of these pa- 
rameters and validate the model. The primary objective 
was to reproduce the known thermal power of the RTGs 
by choosing suitable temperature boundary conditions 
that, in turn, had to agree with telemetered temperature 
readings. In the final results, the modeled RTG thermal 
power was within 1% of the known value, while modeled 
RTG fin root temperatures were always within ±2 K of 
the flight telemetry. 

Both spacecraft utilized a louver system for thermal 
management Louvers mounted on bimetallic springs 
opened in response to high interior temperatures, al- 
lowing excess thermal radiation to escape the spacecraft 
more freely. Louver geometry for twelve two-blade and 
two three-blade louver assemblies was integrated into the 
model. Movements of the modeled blade positions (and 
the resulting calculated louver effective emittance) were 
based on the average temperature of the two nodes on the 
edge of the element that corresponded to the physical lo- 
cation of the actuator housings for each louver assembly. 

In addition to the internally generated heat, the Sun 
was also a significant source of heating, particularly at 
the smaller heliocentric distances. Since the spacecraft 
was facing the Sun, the solar energy was absorbed pri- 
marily by the HGA which largely shadowed the rest of 
the spacecraft from direct solar irradiation except for 
the RTGs. The solar effect became evident as the ab- 
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TABLE II: Error budget for the Pioneer 10 thermal model; 
the contributions shown in percentages, relative to ap. 



Description 



Error 



TABLE III: Doppler residuals (in mHz) after incorporating 
a thermal model, varying the efficiencies of RTG (J7rtg) and 
electrical (rjcicc) heat conversion into a recoil force. 



\y Lllcl oU Lli l^cb 

Quantization error in telemetry data 


^ n 1 % 
2.2% 


^clcc 


0.0000 


Vrtg 

0.0096 0.0120 0.0132 


0.0144 


0.0156 


0.0288 


11 1-W radioisotope heater units 


2.2% 


0.000 














5.57 


Discrepancy in TWT telemetry 


2.5% 


0.295 




10.80 


9.03 




7.37 






Inaccuracy of geometric model" 


5.0% 


0.369 




8.75 


7.11 




5.70 






Modeling of instrument openings' 


5.0% 


0.443 




6.87 


5.50 




4.62 






Subtotal: 


8.1% 


0.480 








4.57 


4.45 


4.54 




RTG surface degradation 


25.0% 


0.517 










4.57 






TOTAL: 


26.3% 


0.886 


4.52 















"Assuming a fore-aft positional uncertainty of 2.5 cm. 

''Heat escaping through instrument openings is ill-documented. 



sorbed energy increased the HGA temperature and was 
then emitted as IR radiation. 

The model was evaluated for Pioneer 10 at eleven dif- 
ferent heliocentric distances ranging from 3 AU to 80 AU 
(see Fig. Four temperature telemetry locations on 
the equipment compartment and two on the instrument 
compartment were examined. The deviation between the 
predicted values from the model and the flight telemetry 
was tracked in each of the eleven cases. The differences 
between telemetered temperature values at the six plat- 
form temperature sensors on board vs. values computed 
from the model reveal a root-mean-square (RMS) error of 
5.1%. This discrepancy is consistent with, and is likely 
a consequence of, the known sources of error listed in 
Table |lll Of these, the most significant is the effect of 
possible surface degradation of the sun-facing RTG sur- 
faces. 

The RTGs were coated with "three mils of zirconia 
(Zr02) in a sodium silicate binder" 



12|. No information 



is available in the literature about the performance of this 
particular type of paint when exposed to solar radiation, 
especially at the relatively high temperatures present on 
the RTG outer surfaces. Similar paints [l3| have experi- 
enced both an increase and a decrease of up to 5% in in- 
frared emissivity. Approximately 25% of the RTG coated 
surfaces were exposed to solar irradiation. A calculation 
that takes into account the relative contribution of RTG 
heat to the total anisotropy yields a corresponding error 
figure of 25% in the overall error budget. 

Other, significant error sources (see Table lll| include 
the FE modeling uncertainty (obtained by comparing the 
FE model output to redundant telemetry), insufficient 
detail in the documentation about the amount of heat 
escaping through instrument openings, and limited accu- 
racy in the available engineering drawings. 

The thermal power Q emitted in a certain direction is 
related to the acceleration ath resulting from the corre- 
sponding thermal recoil force as ath — V {Q/mc), where 
m is the mass of the object being accelerated and c is the 
velocity of light. The dimensionless quantity rj is an "ef- 



ficiency factor" ; its value is 1 if the beam is collimated in 
the direction of the acceleration. It was recognized ear- 
lier ll| that the two main contributions to the thermal 
recoil force are the waste heat generated by the RTGs, 
Qrtg, and the waste heat produced by electrical equip- 
ment inside the spacecraft compartments, Qcicc, yielding 
the anomalous acceleration 



ap = r]rtg{Qitg/mc) + 77clcc(Qclec/"lc) 



(1) 



Our earlier investigation 11] found that the correspond- 
ing efficiency factors, ?7itg and r/dcc, depend only on the 
geometry of the spacecraft, and remain constant over 
time; this result was confirmed by the present analysis. 
The thermal recoil force is well fitted by the model equa- 
tion using ?7itg = 0.0104 and rycioc = 0.406, with a RMS 
error of only 0.78 W, much smaller than the inherent 
quantization error present in the telemetry. 



DOPPLER ANALYSIS AND COMPARISON 

The parameterized thermal recoil force model pre- 
sented in Eq. ([1]) can be incorporated easily into the 
force model used for Doppler analysis. This presented 
us with the opportunity to go beyond simply using a re- 
coil force estimate in trajectory modeling; the Doppler 
analysis can also be used to estimate the values of the 
parameters ryrtg and rydcc completely independently from 
the thermal analysis. 

Results of this part of our analysis are shown in Ta- 
ble mil The spacecraft trajectory was calculated using a 
variety of values for r^itg and rjdcc- The Doppler anal- 
ysis yielded the lowest residuals at ryitg = 0.0144 and 
??cicc = 0.480; the residuals, parameterized by rj^tg and 
770I0C, lie on an elongated elliptical paraboloid surface. 
Consequently, the Doppler analysis is sensitive to the 
overall magnitude of the recoil force: a ~ 20% change in 
the overall magnitude nearly doubles the residual. How- 
ever, the Doppler analysis cannot disambiguate between 
the RTG and electrical contributions to the recoil force: 
an increase in one accompanied by a decrease in the other 
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FIG. 3: Comparison of the thermally-induced and anomalous 
accelerations for Pioneer 10. The estimated thermal acceler- 
ation is shown with error bars. The stochastic acceleration 
estimate from ^ appears as a step function. For reference, 
the Doppler residuals of the stochastic acceleration are also 
shown in the bottom panel. Inner and outer error bars corre- 
spond to the subtotal and total shown in Table HIl 

results in a small or insignificant increase in the Doppler 
residual, as can be seen from the off-diagonal elements in 
Table nil 

The uncertainty of the thermal recoil force estimate 
is large. The magnitude of the estimated thermal re- 
coil force is smaller than the observed acceleration, but 
the stochastic acceleration estimate derived from Doppler 
data is within la (see Fig. An error ellipse of the 
thermal recoil force estimate in the r^rtg — Vcicc parameter 
space of Eq. ([T]) is shown in Fig. |4l This error ellipse 
takes into account the fact that the largest source of er- 
ror (RTG coating) affects only i]^tg- Fig. |4] also shows the 
error of the stochastic acceleration, estimated using the 
best-fit Doppler residual as the la noise floor. The ori- 
entation of this error ellipse corresponds to Table Hill and 
reflects the fact that the using the Doppler data alone, 
the RTG and electrical contributions to the thermal re- 
coil force cannot be well distinguished. 

Lastly, we mention that both the thermal recoil force 
and the Doppler data can be well modeled using an ex- 
ponential decay model in the form a = [2~(*^*''^/'^] cq. 
Using to = January 1, 1980, the best fit parameters 
for the Doppler data are [1] r = (28.8 ± 2.0) yr, oq = 
(lO.lil.O) X lO^^*' m/s^. In contrast, the calculated ther- 
mal recoil force can be modeled, with an RMS error of 
0.1 X 10"^" m/s^, using the parameters r — 36.9 ±6.7 yr, 
ao = (7.4 ±2.5) X 10-1" m/s2. 



CONCLUSIONS 

We presented results from a thermal analysis of the 
Pioneer 10 spacecraft and an independently performed 
analysis of their trajectories using Doppler radio-metric 
data that also incorporated a parameterized on-board 
force model. 
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FIG. 4: Comparison of the RTG (r/rtg) and electrical (rydcc) 
efSciency factors estimated independently from the Doppler 
and thermal analysis using Eq. ([T]). The la error ellipse for 
the thermal analysis corresponds to the values in Table [III For 
the Doppler analysis, the error ellipse corresponds to values of 
77rtg and rycioc in Table Hill for which the excess in the Doppler 
residual exceeds the best-fit residual. 

The comprehensive thermal model yields a recoil force 
characterized by the parameters yy^tg — 0.0104 and 770I0C — 
0.406. The Doppler analysis yielded the lowest residual 
at 77i-tg = 0.0144 and 770100 = 0.480. Numerically, the ther- 
mal analysis based estimate of the recoil force is ~80% 
(of which -35% is produced by the RTGs and -45% by 
electrical heat) of the magnitude estimated from Doppler 
analysis. 

To determine if the remaining 20% represents a sta- 
tistically significant acceleration anomaly not accounted 
for by conventional forces, we analyzed the various error 
sources that contribute to the uncertainties in the accel- 
eration estimates using radio-metric Doppler and ther- 
mal models. On the Doppler side, a significant noise 
fioor exists which is believed to be due to systematics (in- 
cluding the interplanetary charged particle environment, 
etc.) For the thermal analysis, the biggest source of un- 
certainty is the unknown change in the properties of the 
RTG coating, which results in differences in fore-aft emis- 
sivity and an additional contribution to the recoil force 
in the spin axis direction. When we plot these uncertain- 
ties in the T^rtg, '7oioc parameter space, the la error ellipses 
overlap. Wc therefore conclude that at the present level 
of our knowledge of the Pioneer 10 spacecraft and its tra- 
jectory, no statistically significant acceleration anomaly 
exists. 

In closing, we must briefly mention additional avenues 
that may be explored in future studies. First, the case 
of Pioneer 11 was not analyzed at the same level of de- 
tail, albeit we note that spot analysis revealed no sur- 
prises for this spacecraft. Second, the question of the 
anomalous spin-down of both spacecraft remains unad- 
dressed, even though it is plausible that the spin-down 
is due to heat that is reflected asymmetrically off instru- 
ment sunshades. Third, Fig.l2lis strongly suggestive that 
the previously reported "onset" of the Pioneer anomaly 
may in fact be a simple result of mismodeling of the so- 
lar thermal contribution; this question may be resolved 
with further analysis of early trajectory data. Fourth, an 



5 



as-yet unaddressed issue is the possibility of outgassing 
from surface materials, which was shown to have a poten- 
tially observable contribution to the anomalous accelera- 
tion [3l • Fifth, our understanding of systematics in the 
Doppler tracking data can be improved by a detailed au- 
tocorrelation analysis. Sixth, the properties of the RTG 
paint are, in principle, measurable by a thermal vacuum 
chamber test of a hot RTG analogue. Finally, yet an- 
other redundant data set exists in the form of Deep Space 
Network signal strength measurements, which could be 
used to improve our understanding of the spacecraft's 
precise orientation. Nonetheless, it is unlikely that any 
re- analysis of the Pioneer 10 and 11 data set will alter 
our main conclusion: the anomalous acceleration of these 
spacecraft is consistent with known physics. 
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